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We have measured the Pancharatnam relative phase for spin-1/2 states. In a neutron polarimetry 
experiment the minima and maxima of intensity modulations, giving the Pancharatnam phase, were 
determined. We also considered general SU(2) evolution for mixed states. The results are in good 
agreement with theory. 
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I. INTRODUCTION 



In recent years much attention has been paid to the 
concept of geometric phase. Since its discovery by Berry 
p] the subject was widely expanded and subdued to sev- 
eral generalizations J2J, e.g. non adiabatic and non- 
cyclic evolutions as well as the off-diagonal case 0. 
Ever since, a great variety of experimental confirmations 
of the Berry phase and its peculiar properties have been 
accomplished (see e.g. HASHES Ell)- In 1956 Pan " 
charatnam defined <fr = arg(^olV') as the phase acquired 
during an entirely arbitrary evolution of a wave-function 
[l2|. Provided the evolution takes place under condi- 
tion of parallel transport, <p can be identified with the 
noncyclic geometric phase. Quite recently, also a con- 
cept of mixed state phase was developed by Sjoqvist et 
al. [Tflj due to its importance in quantum computation. 
The theoretical predictions have been tested by Du et 
al. 01 an d Ericsson et al. ^3 using NMR and single- 
photon interferometry, respectively. In this letter, we 
report a neutron polarimetry experiment for measuring 
the Pancharatnam relative phase of a spinor subdued to 
an arbitrary SU(2) evolution, implementing the method 
described by Wagh and Rakhecha in 0] . We test its ex- 
tension to the mixed state case put forward by Larsson 
and Sjoqvist [FT! . 



II. THE PURE STATE CASE 
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FIG. 1: Sketch of the neutron polarimetry setup for Pan- 
charatnam phase measurement including evolution of polar- 
ization vector for one specific relative position of spin turners 
at distance Lo and SU(2) coils. 
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In neutron polarimetry the unitary, unimodular operator 
U(a) is used for description of spin rotations (transforma- 
tions) carried out within magnetic fields. By comparison 
of Eq.Q to the matrix representation of U(cE), with 



We consider the experimental setup shown in Fig. ^ A 
monochromatized neutron beam, polarized along the +z 
direction experiences a tt/2 flip about the x-axis, there- 
fore being denoted as a superposition of the orthogonal 
states | ± z) . The beam travels in an overall guide field 
B z of the length L, which is chosen to be parallel to the 
+z direction. Under the influence of B z the polarization 
vector rotates through an angle rj determined by the dis- 
tance L. The two orthogonal components | ± z) of the 
rotated superposition state acquire opposite phase un- 
dergoing the SU(2) transformation Uq, denoted by the 
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one is lead to sets (£, S, Q of SU(2) parameters, a points 
to the direction of the magnetic field constituting the ro- 
tation axis and \a\ = a denotes the rotation angle, a is 
the Pauli vector operator. The resulting beam passes a 
distance L' along the guide field B z , its polarization ro- 
tating through the associated angle 2n7r — r\ (n= 1,2,...). 
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Lo = L+L' is set in such a way that the accumulated ro- 
tation angle within B 2 equals an integer multiple of 2w, 
Lo = nirhv/\fiB z \. Finally, a — ir/2 flip around the cc-axis 
is applied and the output intensity, calculated as 

I = cos 2 £ cos 2 S + sin 2 £ sin 2 (C + r]) (3) 

for the pure state case, is measured. The phase shift 77, 
implemented by variation of the distances L and L' at 
constant Lo causes intensity modulations, yielding the 
extreme values 

Imin = COS 2 £ COS 2 S, (4) 

Lnax = cos 2 £ cos 2 5 + sin 2 £. (5) 

In our case, the Pancharatnam phase <j) is given by sub- 
stituting |^o) by | + z) and \ip) by Uo| + z). The rel- 
ative phase shift between + z) and Uo| + z), cf> — 
arg(+z|U | + z) = 5 + argcos£ for the pure state case, is 
therefore computed from 



c/> = arccos J _ — . (6) 

III. THE MIXED STATE CASE 

More general, a neutron beam passing the setup ex- 
plained above at any arbitrary degree of polarization r 



along the positive z-axis is described by the density op- 
erator 

P = \$ + r&,), (7) 

with the Pauli spin operator a z and < r < 1. The 
measured intensity for this mixed state is given by 

V = i^I+rl, (8) 

which obviously reduces to Eq. © for r = 1 . The relative 
mixed state phase $ 0] is 

$ = arctan[r tan((5 + argcos£)]. (9) 

With the associated extreme values of Eq.© for r > 0, 

Cin = ^+rcos 2 £cos 2 £ (10) 

Cax = iy^+r(cOS 2 ecOS 2 ,5 + sin 2 £), (11) 

one obtains the following expression for the mixed state 
relative phase [l7j : 



* = — V ^/sa+^j+c:- 1/2(1 -.r (12) 



This formula is consistent with the pure state case since 
for r = 1 it reduces to Eq.©. The density operator 
formalism is subject to an ambiguity, in a sense that cer- 
tain inherently distinct physical situations are described 
by identical density matrices. This is often referred to 
as 'decomposition freedom' (see e.g. 18]). Moreover, it 
is not apparent from Eq.Q by what means some mixed 
state of the system should be produced. Considering this, 
we introduced a ir flip in front of the actual setup, result- 
ing in a beam polarized in the — z direction. The inten- 
sities I g and I on , corresponding to 'spin flipper off' and 
'spin flipper on', respectively, were measured for equal 
time intervals. A density matrix was calculated from a 
weighted sum of I ff and I on referring to a certain degree 
of polarization r > and is therefore equivalent to the 
associated mixed state. Since a real polarizer does not 
work perfectly, it is important to note that the incident 
neutron beam can already be considered to be in a mixed 
state, as it was done in this work. For many experiments, 
in contrast, it is sufficient to assume the incident beam to 



be purely polarized in some arbitrarily chosen direction. 

IV. EXPERIMENT 

The experiment was carried out at the tangential beam 
port of the 250 kW TRIGA research reactor of the 
Atomic Institute of the Austrian Universities, Vienna. 
The neutron beam of mean wavelength A = 1.99 A , 
incident from a pyrolytic graphite monochromator, was 
polarized along the +z direction by reflection from a bent 
Co-Ti supermirror array. All spin rotations were imple- 
mented by Larmor precessions about the magnetic field 
axes of DC coils made of anodized aluminium wire wound 
on frames with rectangular profile. After passing the DC 
7r flipper the beam encounters the first tt/2 spin turner 
creating a superposition of the orthogonal states | ± z). 
The polarization vector precesses in the xy-plane through 
the angle r\ around the +z direction in the guide field B 2 
and length L at an angular frequency to — ^fB 2 . The 
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FIG. 2: Rotations on the Poincare sphere associated to parameter sets A: 8^- A \ C^ A ') — > (1.71,0.38,-1.46) rad and B: 

(£ (B) ,<5 (B) ,C (B) ) (1.06,0.17,-1.40) rad for one specific relative position of spin turners and SU(2) coils. 



guide field of length Lo was realized by two rectangular 
coils of 150 cm length along the beam trajectory and a 
width of 12 cm. They were arranged in Helmholtz ge- 
ometry to provide for optimal homogeneity. The com- 
bination of the guide field B z measured by a Hall probe 
(= 5.893 ± 0.022 G) and the distance 4 L (about 46 
cm) was chosen in such a way that more than one pe- 
riod (ij > 2w) of intensity oscillations could be scanned 
for its extreme values. The subsequent SU(2) transfor- 
mation Urj(£, 6, C) was carried out by two coils with mu- 
tually orthogonal orientation wound on the same frame. 
They define an arbitrary magnetic field axis in the xz- 
plane, whose effect on the spin state is denoted by Eq.(J2J. 
Along the distance L' from the SU(2) coils to the second 
spin turner the polarization vector precesses by an an- 
gle 2n7r — r\ about the guide field B z before being turned 
through — 7r/2 about the x-axis, analyzed in the +z di- 
rection by a second supermirror array and finally counted 
by a 3 He detector. The variation of the phase shift r\ is 
implemented by physical translation of the two 7r/2 spin 
turners at constant distance Lo = L + L\ This can be 
done conveniently by steps of some mm, exhibiting the 
intensity oscillations evoked by the phase shift r\. 

Here, we present the results for two settings of SU(2) 
coil currents. For each set of transformation param- 
eters (£, S, £), three minima and maxima of I ff and 
Ion were measured in one scan. Given uncertainties 
contain both systematic and statistical errors. The 
first example is a rotation as depicted schematically 
in Fig. The associated SU(2) parameter set is 

(£( A ),<$( A ),C (A) ) -> (1.71,0.38,-1.46) rad. By applica- 
tion of Eq.@ one obtains the theoretical result <E>{ A ' ) = 
0.37 rad, whereas the value computed from the mea- 
sured data using Ea.fEty is <I> r (A) = 0.38 ± 0.08 rad at 
initial polarization Tq — 0.976 ± 0.004. The measured 
intensities I g, I on and the calculated mixed state inten- 
sities for three different degrees of polarization are shown 
in Fig. [3^ and b, respectively. The second example 
is depicted in Fig. [2J3, being described by parameters 
(f (B >,M B \C (B) ) -> (1-06,0.17,-1.40) rad. For this set 
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FIG. 3: a) Measured intensities I fj and Ion for the parameter 
set (£ (A) ,<5 (A) ,< (A) ) -> (1.71,0.38,-1.46) rad at r ~ 0.98. 
b) Calculated mixed state intensities for three degrees of po- 
larization r. Second order fits (see text) are included. 



(B) — 

the theoretical prediction is = 0.17 rad. The exper- 



imental result is $ 



:b) 



0.16 ± 0.06 rad at initial polar- 



(B) 

izationr ' = 0.981 ±0.005. Measured and calculated in- 
tensities are shown in Fig. 0Jl and b. For the density ma- 
trices derived from the experimental data the results are 



(A) 



rad, $ 



(A) 
0.3 



0.35 ± 0.06 (0.31) rad, $ A 6 j = 0.21 ± 0.06 (0.24) 
0.05 ± 0.07 (0.12) rad for the transforma- 
tion (£ (A) ,<S (A) ,C (A) ) and <4 B s = 0.12 ±0.06 (0.14) rad, 



^ = 0.10 ± 0.07 (0.10) rad, $^3 = 0.02 ± 0.08 (0.05) 
rad for {£}- B \ <5^ B -', C^). The values in parenthesis are 
the corresponding theoretical predictions for $ calculated 
from Eq. © . 
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FIG. 4: a) Measured intensities I g and I on for the parameter 

set (£( B ),(5( B \C (B) ) -» (1.06,0.17,-1.40) rad at r ~ 0.98. 
b) Calculated mixed state intensities for three degrees of po- 
larization r. Second order fits (see text) are included. 

V. DISCUSSION 

The extreme values of the intensity modulations show 
slight irregularities in height as can be seen in the graphs. 
This can be explained as an effect of second order neu- 
trons (wavelength A/2). The monochromatization of the 
beam is governed by the Bragg condition 2dsin#B = nA, 
that - for a certain angle 9b - is fulfilled for n= 1, 2, ... . 
Consequently, the beam used for the experiment is con- 
taminated by a percentage of neutrons with only half the 
intended wavelength. In a time of flight (TOF) measure- 
ment, this percentage was found to be < 7.2 % of the 
overall intensity. Neutrons that travel at double veloc- 
ity spend only half the time in magnetic fields along the 
beam trajectory. Therefore, their spin rotates through 
only half the originally intended angles, showing a be- 
havior rather different from the first order intensity. To 
compute the relative phase, values for I max and I m i n were 
determined from a least-squares fitting model (also shown 
in Fig. |2| and Fig. 0}, taking into account the second 
order neutrons. The careful reader will point out that 
randomness is missing in the mixing procedure described 
above: the relative frequencies for finding the system in 
either one of the pure states \+z) or |— z) are not only well 
known, but even determined by the experimenter. As a 



matter of fact, the actual system, the neutron beam, is at 
no time of the experiment in one of the mixed states as- 
sociated with r=0.8, 0.6 or 0.3. We are aware of the fact 
that the procedure carried out does by no means consti- 
tute a general mixing method, but represents a special 
approach within the wide scope of possible and accepted 
techniques for mixed state preparation. Nevertheless, the 
density matrices are generated by using the measured 
data resulting from the conducted experiments, fulfilling 
the predictions developed by Larsson and Sjoqvist |17j . 
A neutron optical experiment showing the r dependence 
of the phase by application of some other preparation 
method will be carried out in the future. 

The obtained experimental results are in good agree- 
ment with theoretical predictions for the relative mixed 
state phase, although systematic errors occur due to in- 
herent difficulties of the experiment. Even though the 
guide field was constructed carefully, some inhomogene- 
ity (~ 0.4%) along its length cannot be avoided com- 
pletely, accumulating an error of ~ 5° in the angle i] 
along the distance 4 Lo. Furthermore, it was difficult 
to compute values for the readjustment of 4 Lo after a 
change of SU(2) coil currents. Consequently, the cor- 
rect distance between the spin turners was tuned by its 
variation at minimum intensity relative position of spin 
turners and SU(2) coils, thereby exhibiting the maximum 
of visibility. Another point is that the two 7r/2 spin turn- 
ers were situated within the guide field. Therefore they 
did not exactly carry out the designated rotations of the 
polarization vector to the -y direction and back to the +z 
direction. These inaccuracies are intrinsic to this exper- 
iment and can hardly be avoided at considerable effort. 
Further discussions gave rise to interesting ideas for pos- 
sible improvements of the experimental setup concerning 
accuracy and robustness. At present, they are tested and 
will be described in a forthcoming paper |l9j |. 



VI. CONCLUSION 

To summarize, we have measured the Pancharatnam 
relative phase for the mixed state case in a neutron po- 
larimetry experiment. The dependence of the phase of 
the degree of polarization was indicated by calculation 
of weighted sums of the density matrices corresponding 
to the measured intensities. The experimental data is in 
good agreement with theory. Deviations of the predicted 
behavior were explained to arise due to contamination of 
the beam by second order neutrons. 
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